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ABSTRACT 

Recent work by several groups has established the properties of the dwarf satellites to 
M31. We reexamine the reported kinematics of this group employing a fresh technique 
we have developed previously. By calculating the distribution of a \ statistic (which 
we define in the paper) for the M31 system, we conclude that the total mass (disk plus 
halo) of the primary is unlikely to be as great as that of our own Milky Way. In fact 
the x distribution for M31 indicates that, like NGC 3992, it does not have a massive 
halo. In contrast, the analysis of the satellites of NGC 1961 and NGC 5084 provides 
strong evidence for massive halos surrounding both spiral galaxies. 
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'There are reasons, increasing in number and quality, to be- 
lieve that the masses of ordinary galaxies may have been 
underestimated by a factor of 10 or more' (O striker, Peebles 
and Yahil, 1974). 

'...there is no independent observational dynamical evidence 
to support the view that the masses of galaxies are very much 
larger than those derived in investigations of single galaxies 
which have been carried out over many years' (Burbidge, 
1975). 



1 INTRODUCTION 

According to Mateo (1998) more dwarf galaxies have been 
discovered as Local Group members since 1971 than during 
the preceding 222 years. A large fraction of these systems are 
bound to M31. Recently, Evans and Wilkinson (2000, EW), 
Evans et al. (2000) and Cote et al. (2000) have independently 
estimated the total mass associated with the M31 (disk plus 
halo) and found that it is probably less massive than our 
Milky Way. Here we present an alternative way of estimating 
the total mass of M31 using these same data but employing 
a novel technique we have developed (Erickson, Gottesman 
and Hunter 1987, EGH1; and 1999, EGH2). 

The use of the kinematics of dwarf satellite galaxies 
to investigate the masses and extents of galaxy halos is a 
very powerful technique (see EGH1, EGH2, Zaritsky et al. 
1993, ZSFW1 and 1997, ZSFW2). Almost all studies have 
been limited to describing the average properties of a generic 
galaxy because very few individual systems posses enough 
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satellites to allow unique inferences. Owing to the limita- 
tions imposed by projection effects, comparisons between 
observations and models (see EGH2, and §3.1 below) in- 
dicate that about seven satellites are required before the 
signature of a massive halo is likely to be observed. We will 
discuss this further below. In EGH2 we studied the statistics 
of over 71 satellite-primary pairs. Unfortunately, very few 
of the primaries had as many as three satellites. However, 
NGC 1961 (Gottesman, Hunter and Shostak, 1983; hereafter 
GHS) and NGC 5084 (Carignan et al, 1997), as individual 
galaxies, meet the requirements for a substantial number of 
satellites (at least eight for NGC 1961 and seven for NGC 
5084). In the Local Group, M31 is surrounded by a cloud 
of dwarf systems (see EW, for example) and it is our pur- 
pose to examine the kinematics of its satellites using the 
methods we have developed in our earlier papers, and to 
compare the properties of M31 with those of NGC 1961 and 
NGC 5084. In the sections that follow, we will briefly re- 
view our methodology, discuss the data for the three galax- 
ies being compared and draw some conclusions about the 
properties of M31. 



2 METHODS 

We have extended the method introduced by van Moorsel 
(1982, 1987) for the study of binary galaxies to the study 
of massive primaries with dwarf satellites in orbit (greater 
detail can be found in EGH2). We treat each group as if it 
were a set of multiple satellite-primary pairs. We estimate 
the large scale mass, M, of the spiral (obtained from the 
orbital properties of each satellite, which we will call the 
orbital mass) in terms of the dimensionless function \ as 
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M x = r p V r 2 /G (1) 

where r p and V r are, respectively, the projected radial sep- 
aration, and the difference between the projected radial ve- 
locity of each satellite and primary. For r p , we use the chord 
length distance that separates the satellite from M3f rather 
than the arc length. Thus, r p — 2 Rsin(<j>/2). R is the dis- 
tance to the satellite and <j> is the angle between M3f and 
each satellite calculated from standard spherical trigonom- 
etry (Smart, 1977). The data for NGC 1961 and 5084 in 
this paper were obtained from a literature search using 
the NASA/IPAC Extragalactic Database (NED) and differ 
slightly from EGH2. 

If each object is idealized as a point mass, bound to 
the primary (as discussed in EGH2), the function \ (which 
measures the magnitude of the projection effects) depends 
on the true anomaly, v, the angle between the line of nodes 
and the major axis of the orbit, w, the inclination of orbital 
plane, i, and the eccentricity of the orbit, e, 

. 2 ,[cos(u + uj) + ecosu] 2 . 2 . 2 i 

X = sin i II — sin (y + ui)sin i\ 2 . (2) 

1 + e cosv 

As an estimate of the importance of any halo we have 
the observed value, Xobs, 

Xobs = r p V? /G(mi + ra 2 ) (3) 

where mi and m 2 are the masses determined from the rota- 
tion of the individual primary and each satellite respectively. 
The sum of these masses is an estimate of the total mass on 
a small scale. For our systems, m 2 /mi -C 1; hence, m 2 can 
be ignored and the subscripts dropped. Xobs is related to x 
by 

Xobs = {M/m) X - (4) 

Thus, if the mass determined from the rotation curve, 
m (which we will call the rotation mass) , is the total mass of 
the galaxy, the two estimates of this mass will be the same 
and the distribution of Xobs will be the same as that of x- 

If values for the orbital elements in equation (2) are 
chosen randomly, the theoretical distribution of x, P{x)i IS 
found to be sharply peaked for small values of x- That is, 
there is a strong tendency to underestimate the orbital mass, 
M, by almost an order of magnitude, because small values 
of r p and V r are highly favored. Furthermore, as long as the 
rotation curve reflects all the mass of the galaxy, and the 
satellites are bound, P(x) —> as x ~~ *• (1 + e )> which must 
be less than 2. This analysis employs point masses, but in 
EGH2 we showed that this assumption is not critical to the 
conclusion. 

The method has several advantages and we comment 
on two (see EGH2 for details). The ratio, Xobs, is indepen- 
dent of the distance scale, and it compares directly the to- 
tal mass of the primary estimated simultaneously at several 
different radii. Thus we can separate the notion of a massive 
galaxy from that of a galaxy with a massive halo. It is quite 
certain that the spiral galaxies NGC 1961 and NGC 5084 
are very much more massive than M31, but that in itself 
implies nothing about the existence of heavy halos for any 
or all of the systems. We explore this further in the sections 
that follow. 



Table 1. New Primary-Satellite Systems (Ellington, 2001) 



Primary 


Satellites 


r p 


m 


M 


Xobs 


NGC 2532 


CG 0198 


87 


16 


2.9 


0.18 


UGC 1448 


SBS 0155+018 


43 


13 


58 


4.36 


NGC 3677 


CGCG 242-036 


94 


92 


18 


0.2 


NGC 1097 


NGC 1097a 


16 


57 


1.7 


0.03 




ESQ 416-G032 


223 




3.2 


0.05 



Col. 1: name of primary galaxy; Col. 2: name of associated 
satellite(s); Col. 3: projected primary-satellite separation, in kpc, 
based on a Hubble constant H = 67 km s _1 Mpc~ 1 ; Col. 4: 
rotation based mass of primary from HI profile (see EGH2), in 
units of 10 10 Mq \ Col. 5: mass of primary from orbital properties 
of satellites, in units of 10 Mq; Col. 6: the value of Xobs 
■ for each primary-satellite pair. 



Table 2. Primary Galaxies 



Galaxy 


D 




V 7 {^Rmax ) 


Rmax 


Mass (m) 


NGC 224 


0.77 


-301 


235 


30 


38 


NGC 1961 


58 


3934 


409 


36 


140 


NGC 5084 


30 


1726 


328 


68 


170 



Col. 1: name of primary galaxy; Col. 2: with the exception 
of NGC 224 (M31), distance (Mpc) using heliocentric, systemic 
velocity (in units of km s _1 ) (Col. 3) corrected to CMB frame and 
H = 67 km s _1 Mpc -1 ; Col. 4: rotation velocity measured at 
maximum radius (in kpc) observed (Col. 5); Col. 6: the resultant 
mass using a Keplerian calculation, in units of 10 10 Mq. 

3 THE GALAXIES 

3.1 NGC 1961 and NGC 5084 

In Figure 1 we show the observed distribution of Xobs, 
N(x), for the combined, independent, published sets of 
data (EGH2 and ZSFW2), augmented by five new satellite- 
primary pairs found by Ellington (2001). In all cases the pro- 
jected separations of the satellites are less than 300 kpc (we 
adopt H — 67 km s _1 Mpc -1 ). The properties of Elling- 
ton's groups are shown in Table 1. A small fraction, 14 per 
cent, of the 71 data points have Xobs > 2. If contamination 
by non-physical pairs is not a problem and all the satellites 
are bound, then it is these few points (Xobs > 2) that pro- 
vide information about mass distributions that extend well 
beyond the radius of the disks (halos). 

In EGH2 we made extensive comparisons of our data 
with those of ZSFW1 and ZSFW2 (ZSFW1 is really a sub- 
set of ZSFW2 and we will ignore it in this discussion) for 
which the data were chosen using criteria very different from 
those used in EGH1 and its successors. The question we ask 
is: what is the probability of finding a galaxy-satellite pair 
with Xobs greater than two? In EGH2, from the analysis of 
independent data sets, we found that between 13 per cent 
and 21 per cent of the satellites had Xobs values in excess of 
two. We argued that these differences were not systematic 
but were measures of precision. Furthermore, as we discuss 
below in §4.1, the data that followed from EGH1 could be 
modelled marginally by a generic spiral galaxy with halo. 
The details of the model are given in §4.3 of EGH2. A param- 
eter of interest to our discussion here is the eccentricities of 



© 0000 RAS, MNRAS 000, 000-000 



On the Mass of M31 3 




Figure 1. The N(\) distribution for 71 primary satellite pairs discussed in EGH2 augmented by Ellington (2001). The ordinate is N(\); 
the abscissa is Xobs m units of 0.25. The data are restricted to r p < 300 kpc. Ten points show Xobs > 2. One point is off the scale with 
Xobs = 10.6. 



the satellite orbits which ranged from 0.5 to 0.9. The model, 
appropriately defined, was just successful in reproducing the 
observations, but qualitatively the initial peak was low and 
the tail was short. Thus the 21 per cent of the values we 
observed with Xobs > 2 is probably at the high extreme and 
the 14 per cent we found for the combined data we show in 
Figure 1 is a conservative but not unrealistic estimate. On 
this basis, we assume a probability of one in seven. Thus, 
primary galaxies must have about seven satellites before we 
can have a reasonable expectation of finding one with a large 
value of Xobs ■ A probability as large as one in five would cre- 
ate an even greater problem in consideration of M31, as we 
will show. (Note, in EGH2 we discussed a probability distri- 
bution, here we discuss counts, N(x)-) 

Two galaxies that we have studied that meet this crite- 
rion of having at least seven satellites are NGC 1961 (GHS, 
1983) and NGC 5084 (Carignan et al. 2000). As we have 
noted, the essence of our method requires two independent 
methods for estimating the total mass of a galaxy. (1) We 
must be able to measure the rotation velocities of the disk. 
This provides a mass at relatively small radii. (2) We must 
be able to measure the systemic velocities of the satellites 
to the galaxies in question. This enables us to estimate the 
halo masses at larger radii. 

For NGC 5084, the properties of the disk were measured 
by Gottesman and Hawarden (1986) and the satellites were 
investigated by Carignan et al. (1997). For NGC 1961, the 
HI was extensively observed by Shostak et al. (1982). How- 
ever, the nature of the major axis data made it difficult to 
estimate V(Rmax)- EGH2 employed a value consistent with 
the lower range of the Dutch data. Here, in addition, wc 
consider the width of the total, or global, HI spectrum. We 



Table 3. NGC 1961 Satellites and Measured Xobs Values 



Galaxy 


AV 


r p 


M 


Xobs 


CGCG 329-011 


+174 


129 


91 


0.65 


CGCG 329-009 


-62 


152 


14 


0.10 


LEDA 138826 


-134 


188 


78 


0.56 


UGC 3342 


+40 


231 


8.6 


0.06 


UGC 3344 


+348 


342 


960 


6.88 


UGC 3349 NOTES01) 


+358 


456 


1357 


9.72 


CGCG 307-021 


+67 


464 


48 


0.35 


UGC 3349 


+380 


522 


1749 


12.54 


CGCG 329-010 


+82 


528 


82 


0.59 



Col. 1: name of the satellite galaxy taken from NED; Col. 2: 
radial velocity difference (satellite-primary), in units of km s -1 ; 
Col. 3: projected separation, in kpc; Col. 4: the resultant orbital 
mass using data in columns 2 and 3, in units of 10 10 Mq; Col. 
5: Xobs = (M/m)x, the ratio of the mass determined from the 
kinematics of the satellite in col. 4 to the mass of the primary 
from Table 2. 

find that the maximum velocity lies within a range of 367 
and 450 km s _1 . Adopting V(R ma x) = 409 km s _1 yields a 
mass of 1.40 x 10 12 Mq (see table 2). This is greater than 
used in EGH2 and produces smaller values of x- 

In Figures 2 and 3 we show the N(x) distributions for 
these two galaxies; the data for all satellites are plotted. In 
both cases, and more so if the data are combined, there is 
clear evidence for the existence of massive halos with sizes 
greater than 100 kpc in the case of NGC 5084, and greater 
than 300 kpc for NGC 1961. The data for M31 are shown in 
Figure 4. 
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Figure 2. The N(\) distribution for NGC 1961. The data are taken from Tables 2 and 3. The maximum projected separation for the 
satellites is 528 kpc. One point is off the scale with Xobs = 12.54. 
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Figure 3. The N(\) distribution for NGC 5084. The data are taken from Tables 2 and 4. The maximum projected separation for the 
satellites is 162 kpc. 
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Table 4. NGC 5084 Satellites and Measured Xobs Values 



Galaxy 


AV 


r p 


M 


X.obs 


Gl 


-40 


162 


6.0 


0.04 


G2 


-238 


97 


128 


0.75 


G3 


-295 


71 


143 


0.85 


G4 


-635 


94 


877 


5.2 


G5 


+479 


51 


273 


1.6 


G6 


-387 


61 


212 


1.3 


<G7-G8> 


+356 


134 


394 


2.3 


G9 


+167 


103 


67 


0.39 



Col. 1: name of the satellite galaxy taken from Carignan et al. 
(1997); Col. 2: radial velocity difference (satellite-primary), in 
units of km s — 1 ; Col. 3: projected separation, in kpc; Col. 4: 
the resultant orbital mass using data in columns 2 and 3, in units 
of 10 10 Mq; Col. 5: Xobs = (M/m)x, the ratio of the mass 
determined from the kinematics of the satellite in col. 4 to the 
mass of the primary from Table 2. 

3.2 The M31 Group 

The Local Group is dominated by our own Milky Way and 
by M31 (NGC 224). They are separated by 770 kpc and 
each is surrounded by a group of dwarf satellites. Our prime 
references for the properties of M31 and the M31 group have 
been EW, Evans et al. (2000), Cote et al. (2000) and Irwin 
(private communication). However, to ensure uniformity our 
ultimate authority has been NED, which we use throughout 
this paper. (Dr. Wilkinson has informed us that the position 
given by Evans et al. (2000) for And vii is in error; our values 
are correct.) 

In Table 2, we list the properties of the three primary 
galaxies, while in Tables 3 — 5 we list the properties of each of 
the satellite groups. However, before considering the N(x) 
distribution for M31 a few comments are required. 

First, what satellites should be included in the analysis? 
We are concerned in particular with the Pegasus and IC 
1613 dwarfs. These galaxies lie at projected separations in 
excess of 400 kpc from M31, which is more than one-half the 
distance to the Milky Way. We presume the orbital period 
of the satellites to be at most a Hubble time. Employing a 
simple Keplerian calculation, the radius of the largest orbit 
is proportional to M 1/3 . If we adopt the mass of NGC 5084 
(170 x 10 10 Mq) as a fiducial indicator, the maximum radius 
that satisfies this condition is about 350 kpc. Employing 
this size may be in error, but it will allow us to compare 
the three galaxies based on the initial assumption that all, 
including M31, are massive systems. It is unreasonable to 
expect a two body calculation to be applicable for satellites 
at such large distances from the primary (M31) in the nearby 
presence of our own massive galaxy. However, as we will 
show, even if these two distant dwarf galaxies are included 
in the discussion, the M31 group is unusual. 

Second, Andromeda is a highly resolved system with 
an optical diameter in excess of three degrees. This makes 
it unlike any of the other galaxies we have studied. Our 
lines of sight to the centre of the galaxy and to the satel- 
lites are not parallel. Indeed, as independent distances can 
be measured to the satellites, most studies (see for example 
EW) are done in an Andromedacentric system. As noted in 
equation (3), our formalism has been developed for projected 
velocity and position coordinates. We cannot correct for line 
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Table 5. M31 Satellites and Measured Xobs Values 



Galaxy 


A V 

L\ V 


r p 




Xobs 


M32 


+96 


5.4 


1.17 


0.03 


M33 


+75 


198 


25.5 


0.67 


NGC 147 


+118 


100 


32.2 


0.85 


INLr^ loO 


+ 1U t 


yo 


zo.z 


U.OD 


tv t c A /~< one; 


+61 


Q O 

o.z 


A or: 
U.oO 


n no 
U.Uz 


IC 10 


-33 


246 


6.35 


0.17 


*IC 1613 


-60 


520 


43.9 


1.16 


LGS 3 


-30 


266 


5.36 


0.14 


*Peg dwarf 


+85 


412 


68.4 


1.8 


And i 


-85 


44 


7.42 


0.20 


And ii 


+82 


138 


21.8 


0.57 


And iii 


-58 


67 


5.3 


0.12 


And v 


-108 


108 


28.9 


0.76 


And vi 


-65 


264 


25.7 


0.68 


And vii 


+22 


217 


2.53 


0.07 



Col. 1: name of the satellite galaxy taken from NED, IC 1613 
and Peg dwarf (Peg) are marked * as their membership in the 
group is uncertain; Col. 2: radial velocity difference (satellite- 
primary), in units of km s _1 ; Col. 3: projected separation, in 
kpc; Col. 4: the resultant mass using a Keplerian calculation, in 
units of 10 10 Mq; Col. 5: Xobs = (M/m)x, the ratio of the mass 
determined from the kinematics of the satellite in col. 4 to the 
mass of the primary from Table 2. 

of sight effects on the observed radial velocity without know- 
ing something about space motions. We expect the effects 
not to be large and certainly to be less than factors of three 
or more required to alter our discussion. We have compared 
a x'obs calculated using the correct satellite distance from 
M31, as well as Xobs using projected separations and a dis- 
tance to M31 of 770 kpc. For our data, the Xobs values from 
V 2 r p are about 90 per cent of those obtained from V r 2 r. 

To gain some insight on this problem it is instructive 
to consider the special case of a circular orbit of radius r 
and circular speed V. In order to maximize the x values at 
all orbital phases, we let the orbital plane lie perpendicu- 
lar to the plane of the sky (i = 90°). The values of V 2 r p 
range from to V 2 r, with the angle averaged values being 
4/(3tt) V 2 r = 0.4244 V 2 r. It should be noted that the coeffi- 
cient is independent of the ratio r/d, where d is the separa- 
tion between the centres of the Milky Way and M31. A sim- 
ilar calculation of the angle averaged value V?r = 0.5 V 2 r. 
The coefficients 4/3-7T and 1/2 are maximum values because 
7r/2 is the maximum inclination for any orbit. For this sim- 
ple example, the angle average of V r r p is roughly 85 per 
cent of the angle average of V 2 r. 

Although it makes no essential difference to our conclu- 
sions, we will adhere to projected separations in order to be 
as consistent as possible with our earlier studies. In Figure 
5 we show the N(x) distribution for the thirteen satellites 
of M31, for which r v < 350 kpc. This should be compared 
with Figure 6, which shows the combined data for NGC 1961 
and NGC 5084 for which r p < 343 kpc. Unlike these latter 
systems, M31 shows no Xobs values greater than 2. Only if 
we include the observations for the Pegasus dwarf, for which 
r p = 412 kpc, do we get a value as large as 1.8 (see Figure 4). 
Given the projected separation, this is surprisingly small, if 
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Figure 4. The N(x) distribution for M31 as listed in Table 5. Data for M31, the primary, is taken from Table 2. The maximum projected 
separation for the satellites is 520 kpc. 



M31 has a heavy halo. It is considerably less than the Xots 
values of five or more seen for NGC 1961 and NGC 5084. In 
both of these cases we are considering fewer satellites than 
for M31, and, for NGC 5084 in particular, the satellites, on 
average, have smaller projected separations than those ob- 
served for M31. 

The number of galaxies studied by EGH2 was large 
enough that a separation could be made between systems 
having r p < 300 kpc, and those having r v > 300 kpc. The 
first group of 71 points has been shown in Figure 1. For 
the second case, 11 of 37 Xobs values were greater than two 
and of these two had values 10 < Xobs < 21. Based upon 
our studies of bound systems with massive halos, we expect 
at least 14 per cent (one in seven, as noted in §3.1) of the 
satellites to show Xobs values between two and ten, if the 
projected radius of the group is about 300 kpc. None of the 
satellites of M31 that meet this criterion have Xobs ratios 
greater than two. Furthermore, if the radius of interest is 
extended out to 500 kpc there should be larger values of 
Xobs, some approaching a value of 20. However, Xobs = 0.32 
for IC 1613 (r p = 520 kpc). Thus, the data for M31 differ 
significantly from the general galaxies studied in EGH2, and 
of more interest, the system differs substantially from NGC 
1961 and NGC 5084, galaxies which we believe show conclu- 
sive evidence for extended and massive halos. If we consider 
'satellites' orbiting at distances as large as 530 kpc, then 
three out of nine dwarfs associated with NGC 1961 have 
Xobs values greater than two (6.9, 9.7 and 12.5). When we 
include the two points already noted for NGC 5084 (which 
has not been surveyed for satellites with such large orbits), 
a total of five satellites out of a population of 17 indicate 
the presence of a massive halo. None of the satellites for 
M31, out of a population of fifteen, meet this criterion. The 



difference is further emphasized by comparing the rms ve- 
locity dispersions of the satellites to the rotation velocities 
of the primaries. For M31 this ratio is 0.33, while for the 
combined data of NGC 1961 + NGC 5084 the ratio is 0.88 
(for all three galaxies, r p < 412 kpc). There appear to be 
two possible conclusions: (1) there is no massive halo asso- 
ciated with M31, or (2) the characteristics of the orbits of 
the satellites to M31, as a group (including the effects of dy- 
namical friction), are systematically different from those of 
more isolated 'field' galaxies. Either way, the results for the 
M31 system are not consistent with those of other massive 
galaxies. 



4 DISCUSSION 
4.1 Kinematics 

In their discussion on determining the mass of spherical sys- 
tems, Bahcall and Tremaine (BT, 1981) employed the kine- 
matical properties of a collection of bound particles. They 
considered an estimator based on the projected mass 



vlR 
~G~ 



(5) 



where both the velocity and radial separation are projected 
quantities, q serves as an estimator for the mass M of the 
system. For an arbitrary distribution function they showed 
that the expectation value of q becomes 



< q >-- 



32 v 



2 < e 2 >). 



(6) 



If all the particles have circular orbits < e 2 > = 0, while 
for radial orbits < e 2 > = 1. These arguments lead them to 
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Figure 5. The N(\) distribution for M31 as in Figure 4 except that the maximum projected separation for the satellites is 266 kpc. 



define two specialized estimators based on the mean of the 
projected mass < q > — < v^R > /G : 



!6 

Mr = ^< 9 > 



32 

M L = — < q > 

7TG 



(7) 



(8) 



For isotropic orbits < e 2 > = 0.5, and Mi = M. For 
linear orbits < e 2 > = 1, and Ml = M. Equations (6) — (8) 
show in a very simple way, but clearly, that the estimated 
mass of M31 depends on the average eccentricity of the or- 
bits of the satellites. The greater the typical eccentricity, 
the greater the underestimate of the mass of the galaxy by 
< q >, the estimator of the projected mass. From equation 
(6) the variation between extremes can be as much as a fac- 
tor of three. Of course, the analysis above does not consider 
the influence of a halo with a logarithmic potential. 

However, EGH2 did consider the general problem of 
orbits in massive, isothermal halos. As the orbits are not 
closed, the eccentricity was defined, as can be done in the 



,)/(r 



+ r„ 



where 



i ap g and r vrg are the apogalactic and perigalactic radii re- 



Keplerian problem, e 

r a 

spectively. In a pure halo potential the large x values were 
found to originate in the vicinity of perigalacticon from or- 
bits with large r apg (that is, where the apogalactic radius ap- 
proaches that of the halo, r apg — > rjj). In this fashion EGH 
were able to explain the large \ values of about 20 that were 
observed in EGH2 as arising from systems with massive ha- 
los and radii in excess of 300 kpc. The authors found that, 
as the eccentricity of the orbits varied from 0.90 to 0.00, the 
maximum value of \ increased from 9.7 to 24.7. The mass 
of the primary was unchanged for all of these experiments. 
As the orbital eccentricities increase for a given model, the 



estimate of the mass (equivalent to < q >) decreases, as was 
noted by BT. Thus, one way of 'hiding' mass, as indicated 
by the N(x) distribution, would be to have essentially linear 
orbits. While the galaxy halo would be massive, the orbits 
would show an unexpectedly low value for \max- 

In EGH2 we explored a large number of N-body models 
designed to simulate a range of N(x)- In order to compare 
these results with M31 certain adjustments to the models 
are required. The asymptotic maximum velocity has to be 
scaled to 235 km s^ 1 ; the Hubble constant has to be revised 
to H = 67 km s _1 Mpc -1 , and the characteristic eccen- 
tricity of the orbits has to be increased at least to 0.8. (See 
the discussion associated with §4 of EGH2.) If the apogalac- 
tic distance of the satellites is, for example, 350 kpc then a 
larger eccentricity implies that perigalacticon is within the 
disk and the effect of interactions should be severe and clear 
to the observer. This does not appear to be the case in gen- 
eral. (However, M32 is a notable counter example, which we 
will discuss below.) 

The first model we discuss (model 5 in EGH2) was cho- 
sen to reproduce the data of EGH1 and its successors. The 
second model (model 6) was developed to satisfy the data 
of EGH2, but did not reproduce the very largest values of 
Xobs- The third model (model 3) was formulated to describe 
the full range of data in EGH2. 

Model 5 was a one-galaxy model, a generic disk plus 
halo combination. The halo radius was 90 kpc and the total 
mass of the system was 9 x 10 11 M . This model would 
require many of the satellites of M31 to orbit outside the 
galactic halo. 

Model 6 was constructed to model the data in EGH2. 
As we will discuss in more detail below, the full range of data 
in EGH2 required a two-component model, a very massive 
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Figure 6. The combined N(\) distribution for NGC 1961 and NGC 5084. The data are taken from Tables 2, 3 and 4. The maximum 
projected separation for the satellites is 342 kpc. 



disk plus halo constituent, and an essentially haloless disk. 
Model 6 was capable of reproducing the data of EGH2 pro- 
vided the highest values of Xobs were ignored (possible in- 
terlopers) . The massive component of this model had a halo 
radius of 225 kpc, and a total mass of 2 x 10 12 Mq. 

In order for the tail of model N{\) distributions to ap- 
proach the largest values of Xobs found in EGH2, model con- 
stituents with halo radii in excess of 300 kpc were required. 
Model 3 was the most successful example. Two types of 
galaxy were required to model the full range of N(x)- One, 
a disk with a massive halo was required to reproduce the 
high Xobs tail, and two, a disk-only component was neces- 
sary to match the peak in N(x) at small values of Xobs- As 
in model 6 above, we will assume that the massive element 
duplicates the properties of M31. For this model, the halo 
radius was 350 kpc, and the total mass of the system was 
3 x 10 12 Mq. 

Each of these models was developed independently from 
numerical experiments designed to satisfy different criteria. 
However, as the data were common and several of the char- 
acteristics of the models were common (such as the radius 
of any disk), it is not surprising that the mass of successful 
model galaxies scales with halo radius. With regard to M31, 
however, these models have a singular failure. The maximum 
value for x ls m the range of 9 — 12. Such models would de- 
scribe the properties of NGC 1961 and NGC 5084 but not 
M31. In order to approach a maximum value of Xobs — 2, 
the halo masses of the model galaxies must be reduced. 

We conclude, it is very unlikely that the mass of M31 
within a radius of 350 kpc is greater than 6 x 10 11 Mq. This 
solution requires eccentric (near linear) orbits and suggests 
that the halo dynamics of M31 are quite different from those 



of the galaxies studied by EGH2 and certainly are quite 
different from either NGC 1961 or NGC 5084. 

The data for M31 bear a strong similarity to those of 
NGC 3992. We have argued from several perspectives that 
this galaxy is unlikely to have a massive halo. (See EGH2 
page 158 for further references and for a discussion of the 
data.) Like M31, NGC 3992 has a very compressed N(x) 
distribution. 

Courteau and van den Bergh (1999) suggest that IC 
1613 (the other questionable dwarf 'satellite' of M31, lying 
at a slightly larger distance from M31 than Peg) is a free- 
floating member of the Local Group rather than a satellite 
of M31. A similar argument can be made for the Peg dwarf, 
given its large separation from M31 and the nearby pres- 
ence of the Milky Way. If the radius of any putative halo 
for M31 is closer to 250 kpc, the mass is not greater than 
4 x 10 11 Mq. Braun (1991), from HI observations, finds a 
falling rotation curve with a value of 200 km s~ x at a radius 
of 26 kpc (M = 2.4 x 10 11 Mq). No heavy halo is required 
to reproduce his observations. However, if the rotation curve 
remained flat out to a distance of 260 kpc (for which there is 
no observational evidence), the implied mass of an isother- 
mal halo for M31 would be 2.4 x 10 12 Mq. Such a system is 
too massive to explain the nature of the N(x) distribution 
for the satellites. To conclude this part of the discussion, 
we emphasize that it is very difficult (unless nature con- 
spires) to have a galaxy system (M31 + satellites) with a 
total mass larger than 4 x 10 11 Mq, and not have several 
satellites with Xobs > 2 and Xmax substantially greater than 
1.8. If M31 has such a large halo mass, the kinematic prop- 
erties of the satellites should resemble those of NGC 1961 
and NGC 5084, which they do not. 



© 0000 RAS, MNRAS 000, 000-000 



On the Mass of M31 9 



4.2 Dynamical Friction 

EGH2 included dynamical friction in the model calculations, 
and its effects cannot be ignored if the satellite mass ni2 > 
10 9 M Q . In Model 6 of EGH2, N(x) distributions for satel- 
lites having masses of m 2 = 2 x 10 9 M were compared to 
distributions for satellites of m,2 = 0.67 x 10 9 Mq. After 
integrating the orbits of the higher mass satellites for one 
Hubble time, dynamical friction was shown to truncate the 
tail at large \ values, and to increase the peak of the distri- 
bution at low values of \ . In other words, for the galaxies 
with massive halos studied by EGH2, the satellites cannot 
have masses much larger than 10 9 Mq or they will spiral 
into the nucleus of the primary. One uncertain parameter 
controlling deceleration caused by dynamical friction is the 
Coulomb logarithmic coefficient A, which was taken to be 
6.8. One can define an effective satellite mass M sa tA, which 
in our studies was 10 9 Mq , and this value was utilized as the 
mass of the satellites in our numerical experiments. If M31 
has a massive halo the satellites must have small masses, as 
denned above. If they have large masses the existence of a 
massive halo would be very doubtful. Unfortunately, this is 
not a sensitive test as most of the satellite masses (that have 
been measured) are significantly less than 10 9 Mq (Mateo, 
1998), with the notable exceptions of IC 10, M32 and M33. 
The implications for the two groups of satellites need to be 
considered separately. 

For the majority of satellites, the masses are too small 
to be effected by dynamical friction and cannot be impor- 
tant in truncating the high Xobs tail of N(x)- If M31 has a 
substantial halo, the nature of the orbits, and not dynam- 
ical friction, is responsible for diminishing the tail in the 
distribution. 

For the massive satellites, the implications are different 
and interesting. M32 is a special case. It is very close to 
M31 (the radial separation is estimated to be 5 kpc) and 
interactions are important (for example Byrd, 1976; Bekki 
ct al., 2001). At this small separation, this satellite may be 
about to disappear owing to tidal and frictional effects. 

In contrast, both IC 10 and M33 lie at distances from 
M31 greater than 190 kpc. At these distances the frictional 
effect of the halo is unlikely to be large. If the orbits of 
these systems are deeply plunging (as we have strongly ar- 
gued may be the case, if the halo mass of M31 is large) then 
dynamical friction is important and the orbits should have 
decayed over a Hubble time, but they have not. M31 cannot 
have a very massive halo, unless the orbits of these massive 
satellites have a low eccentricity. However, if the orbits of 
all of the satellites are near-circular, then N(x) is not con- 
sistent with the existence of a massive halo. Of course, the 
eccentricity of the satellite orbits could change systemati- 
cally with apogalactic distance, with more circular orbits at 
larger radii. Such complications are not required to explain 
the properties of satellite orbits in the general case, as shown 
by EGH2. Thus, if < e 2 > « 0, the M31 system would be 
most unusual. 

It is clear that the majority of the satellites are not 
massive enough to have suffered significant frictional decel- 
eration regardless of the properties of the halo. The orbits of 
the very few (known) massive satellites should have suffered 
from dynamical friction. If they have eccentric orbits they 
should have decayed and now resemble that of the Large 



Magellanic Cloud (see EGH2) and be on their last few cir- 
cuits before oblivion. On the contrary, they lie at large radial 
distances from M31. Either dynamical friction is insignifi- 
cant or their orbits are not deeply plunging. 

We conclude that: (1) the masses of most of the satel- 
lites are too small to have suffered a significant deceleration 
in a massive halo; (2) the massive and distant satellites may 
not be on highly eccentric orbits if a massive halo exists; 
or (3) most simply, a massive halo does not exist and fric- 
tional decelerations are not expected for any of the satellites 
without regard to location or mass. 

Mass and length are directly proportional to the value 
of H employed in the calculations. The results, taken from 
our models and quoted here, have been scaled to H — 67 
km s _1 Mpc -1 . Our models were developed in EGH2 using 
a Hubble constant of H = 100 km s _1 Mpc -1 . However, we 
investigated the characteristics of these model galaxies for 
different values of H a , and no significant differences in our 
conclusions were found for 50 < H < 100 km s _1 Mpc -1 . 



5 CONCLUSION 

In conclusion, we emphasize that the models for M31, that 
we have considered, were contrived to produce the greatest 
halo mass consistent with the results of EGH2. We find no 
convincing evidence for the existence of a massive halo, nor 
do we require unusual characteristics for the orbits of the 
satellites. 'Essentia non sunt multiplicanda praeter necessi- 
tatum' (see for example Ockham, W. as printed in 1495). 

On the basis of the analysis we have presented, we be- 
lieve that the Keplerian mass, consistent with Braun's 1991 
observations of M31, is close to the total mass of the galaxy. 
Certainly, if NGC 1961 and NGC 5084 are typical of galax- 
ies with massive halos, then M31 is not of this class. The 
N(x) distribution is too compressed. Our limit is consistent 
with, but at the lower bounds of, the models presented by 
EW and Cote et al. (2000). However, we are dealing with 
a small number of satellites and caution is required before 
drawing definitive conclusions. The detection of two or three 
new satellites of M31, with separations less than but veloc- 
ities relative to M31 greater than that of the Peg dwarf, 
might alter our conclusions. In addition, the velocities of 
the satellites must be made more secure. However, in the 
absence of such data we are in complete sympathy with the 
suggestion that the Milky way is the massive member of the 
Local Group (Wilkinson and Evans, 1999). 
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